Flow reversals in turbulent convection via vortex reconnections 
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We employ detailed numerical simulations to probe the mechanism of flow reversals in two- 
dimensional turbulent convection. We show that the reversals occur via vortex reconnection of 
two attracting corner rolls having same sign of vorticity, thus leading to major restructuring of the 
flow. Large fluctuations in heat transport are observed during the reversal due to this flow recon- 
figuration. The flow configurations during the reversals have been analyzed quantitatively using 
large-scale modes. Using these tools, we also show why flow reversals occur for a restricted range of 
Rayleigh and Prandt numbers. 
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growth of the secondary modes at the expense of the 
primary modes is akin to the cessation-led reversals re- 
ported by Brown and Ahlers js] and Mishra et al. fllj j . 
and to the emergence of quadrupolar mode in dynamo 
reversals Q- 

One important question that persists is why do flow 
reversals occur in a restricted parameter regime? In this 
letter, we explain this through a quantitative investiga- 
tion of the convective flow structure for a range of pa- 
rameters using Fourier basis decomposition. This scheme 
allows for accurate representation of the flow [l2| and lets 
us quantify when reversals occur in terms of amplitudes 
of the modes. In addition we show that flow reversals oc- 
cur via vortex reconnections, and connect them to vortex 
dynamics. We investigate in detail the heat transport in 
the system during a reversal and find large fluctuations 
in it as a result of the flow reorganization during the re- 
versals. We restrict our study to 2D flow reversals whose 
flow structures are well represented by experiments con- 
ducted in quasi two-dimensional geometry Q. 

We solve the equations governing Rayleigh-Benard 
convection under Boussinesq approximation 12j in a 2D 
box of aspect ratio 1 with no-slip walls on all sides. The 
side walls are insulating, while the top and the bottom 
walls are maintained at constant temperatures. The 
simulations are performed using the open-source code 
NEK5000 [3] that employs spectral element method. 
We use 28 x 28 spectral element with 7 th order polyno- 
mial, thus we have an effective grid of 196 x 196 points. 
The grid density is higher at the boundaries in order to 
resolve the boundary layers. We perform simulations for 
Pr = 1 and Rayleigh numbers ranging from 10 5 to 10 9 . 

To quantify the flow structures, we project the nondi- 
mensionalized horizontal velocity (u), vertical velocity 
(v), and temperature (T) onto the following basis: 



Several experiments [1H8J and numerical simula- 
tions on turbulent convection exhibit "flow re- 
versals" in which the probes near the lateral walls of 
the container show random reversals (also see review 
articles [HI). These reversals have certain similarities 
with magnetic field reversals in dynamo and Kolmogorov 
flow 0- Researchers typically study convection in a 
controlled setup called "Rayleigh-Benard convection" in 
which a fluid confined between two plates is heated from 
below and cooled at the top. The two nondimensional 
numbers used to characterize the flow are the Rayleigh 
number (Ra) , which is the ratio of the buoyancy term and 
the diffusive term, and the Prandtl number (Pr), which 
is the ratio of the kinematic viscosity and the thermal 
diffusivity. Cioni et al. [l L Niemela et al. Q , Brown and 
Ahlers [3(, and Xi et al. J] performed convection exper- 
iments on mercury, water, and helium gas in a cylin- 
drical geometry and observed reversals for Ra > 10 8 . 
Sugiyama et al. [&] and Vasiliev and Frick [f| studied 
reversals in a rectangular box with water. No reversal 
was observed for a cubical box (aspect ratio 1), but a 
quasi two-dimensional box (aspect ratio < 0.2) exhibits 
reversals for a band of Rayleigh and Prandtl numbers [§[ . 
Surprisingly, a cubical box containing mercury shows re- 
versals [7] , indicating a strong role played by the Prandtl 
number and geometry in the reversal dynamics. 

Several theoretical models have been invoked to ex- 
plain flow reversals. Benzi and Verzicco [9] and Sreeni- 
vasan et al. |14l| used stochastic resonance, while Ara- 
juo et al. [15| employed low-dimensional models with 
noise to explain reversals. Brown and Ahlers |3| and 
Mishra et al. . 1 lj showed that in a cylindrical geome- 
try, the flow reversals are induced by a rotation or cessa- 
tion of large-scale flow structures. For two-dimensional 
box geometry, Sugiyama et al. [8j relate the flow rever- 
sal to the growth of the corner rolls due to the plume 
detachments from the boundary layers. Chandra and 
Verma 12] studied the reversals quantitatively by repre- 
senting flow structures as Fourier modes and showed that 
during the reversals, the amplitude of the first Fourier 
mode (k x — l,k y = 1) becomes very small, while the 
Fourier mode (k x = 2,k y = 2) gains strength. The 
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(v, T) = y^(0 m , Wi T m> „) cos(to7T2;) sm(niry) (2) 



It has been shown that the above basis is a good repre- 
sentative of the flow field in a box geometry 0. We 
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FIG. 1. Steady state flow structures at different Rayleigh numbers: (a) for Ra < 10 5 , a single roll with the dominant 
(ra = 1, n = 1) mode, (b) for 10 6 < Ra < 10 7 , two rolls stacked on top of each other corresponding to the mode (1, 2), and (c) 
for Ra > 2 x 10 7 , two corner rolls with a dominant roll aligned along the 45 degree diagonal, a configuration dominated by the 
modes (1, 1) and (2,2). 



choose the above formalism over proper orthogonal de- 
composition (POD) in order to study the temporal evolu- 
tion of the now structures. The mode with wavenumber 
(m, n) corresponds to a flow structure with m rolls in the 
x direction and n rolls in the y direction. In the follow- 
ing discussion we will use the above basis to analyze the 
mechanism of flow reversals as well as the range of Pr 
and Ra for which reversals take place. 

An important puzzle in 2D turbulent convection is why 
flow reversals take place only for a range of Prandtl and 
Rayleigh numbers. It has been shown that the corner 
rolls (represented by (m = 2, n = 2) mode) are crucial 
for the occurrence of reversals [1, EH. In particular, the 
relative strength of the (2,2) mode with respect to the 
dominant mode determines if a reversal will occur or not. 
We observe three distinct flow structures for the range 
of Ra = (10 5 — 10 9 ) performed in our simulations (see 
Fig. [1]). The mode (1,1), representing the large single 
role, is dominant till Ra ~ 10 5 (Fig. (Ha)). After this, 
the mode (1,2) representing the two horizontally stacked 
rolls dominates till Ra w 10 7 (Fig.QTb)). The mode (2,2), 
which is a dominant player for the flow reversals, is born 
only after Ra « 2 x 10 7 (Fig. [He)), which is the reason 
for the absence of reversals for lower Ra. 

In the third regime, Ra > 2 x 10 7 , the averaged value of 
1^2,2 1 / |^i,i I falls monodically from 0.45 at Ra = 2 x 10 7 to 
0.10atRa = 10 9 (see Fig.EJ). The decrease in |v 2 ,2|/K>i,i| 
is due to the increase of possibly due to inverse 

cascade of energy. This result indicates that for higher 
Ra, the corner rolls become weaker compared to the (1,1) 
mode. Therefore, reversals are observed only for a narrow 
band near Ra = 2 x 10 7 where the (2,2) mode is strong 
enough. We observe similar behavior for Pr = 10 except 
that the (2,2) mode appears for Ra « 10 6 , a value lower 
than that for Pr = 1. Consequently the range of Ra 
exhibiting flow reversals is broader for Pr = 10 compared 
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FIG. 2. Plot of the ratio |«2,2 |/|«i,i I vs - R- a - The mode 1)2,2, 
which quantifies the strength of the corner rolls, is born only 
near Ra = 2 x 10 7 . The average value of the ratio |i>2,2|/|«i,i | 
decreases from 0.45 to 0.10 as Ra increases from 2xl0 7 to 10 9 . 
Reversals are observed for a narrow band near Ra = 2 x 10 7 
where #2,2 is strong. 



that Pr = 1, consistent with the results of Sugiyama et 
al. 0. 

Now we probe in detail the process of flow reversal by 
studying the flow structures of six snapshots during one 
of the reversals (see Fig. A movie of the flow reversal 
can be downloaded from [16[. The starting point of the 
reversal is the stable configuration shown in Fig. QJc). 
At first, the mode (2,2) (or the corner rolls) grows at the 
expense of the mode (1,1), as evident from the subfig- 
ures [3]^a,b). The top-left and bottom-right corner rolls 
have vorticity in the same direction, hence they attract 
each other and come close due to the vortex dynamics 
in 2D [17j. Since the velocities of the streamlines are 
directed in opposite directions, they reconnect, thus con- 
verting two corner rolls into a single roll (see Fig. [3X C )) - 
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FIG. 3. Velocity and temperature profiles for six snapshots during the reversal: (a,b) Growth of corner rolls; (c) The two 
corner rolls at the upper-left and bottom-right corners reconnect to form a large single roll. The streamlines, represented by 
black curves, combine via vortex reconnections; (d,e) The flow reconfigures itself via rotation of the large roll formed after the 
reconnection. There is a strong nonlinear interaction among the modes (1,1), (2,2), (1,3), and (3,1) during these events; (f) 
The flow stabilizes to a quasi-steady state configuration. 



As a result of the reconnection, the new large roll has the 
cumulative vorticity in the same direction. The reconnec- 
tion event leads to a change in the flow topology, and is 
similar to 2D magnetic field reconnections in magnetohy- 
drodynamics 18J. The new large roll aligned along the 
-45 degree diagonal (Fig. [3](c) ) has vorticity opposite to 
that of the large roll before the reversal (Fig. [He)) 

The strength of the new roll grows along with an emer- 
gence of two new secondary modes (1,3) and (3,1), which 
are generated as a result of triad interaction with the con- 
dition k = p + q 0. Note that (1, 3) = (2, 2) + (-1, 1) 
and (3,1) = (2,2) + (1,-1). The subsequent rotation 
of the newly formed roll leads to flow configurations of 
Figs.QJd) and (e), which are dominated by (1,1), (2,2), 
(1,3) and (3,1) modes. Fig.[3Jd) illustrates three horizon- 
tally stacked rolls corresponding to the (1,3) mode. Af- 
ter around 0.01 thermal diffusive time units, the system 
attains a quasi-steady state with the roll aligned along 
—45 degree diagonal with an opposite vorticity (Fig,[3]Jf)) 
compared to the original one (Fig. HJc)). The process 
would repeat for the next reversal with a difference that 
the bottom-left and the top-right corner rolls would re- 
connect during the next reversal. 

The modes (1,1) and (2,2) are the most dominant ones 
for the quasi-steady configurations, with the (1,1) mode 
switching sign between reversals. The sign of the (2,2) 
mode or the sense of rotation of the corner rolls however 
remains unchanged after the reversal, as shown through 
symmetry arguments [l2l |. Another surprising observa- 
tion is that the the amplitude of the (2,2) mode is always 



positive. As a result, the hot plume of the corner rolls (at 
the bottom plate) always ascends via the vertical walls, 
and the cold plume descends via the vertical walls (see 
Figs. [He) and[3Jf)). This seems to be a generic feature 
of many convection simulations and experiments, but its 
mathematical justification eludes us at present. 

We now discuss heat transport during the reversal. 
The global Nusselt number Nu = j[(df(y)/dy) + vT]dr 
exhibits large fluctuations including negative values dur- 
ing a short time interval (see Fig. 0JA)). Here, the 
nondimensionalized T(y) is the averaged temperature 
over the cross-section at height y. The negative Nus- 
selt number occurs for the flow configurations resem- 
bling Fig. OHe) in which the hot fluid parcel in the mid- 
dle of the box descends while the cold parcel ascends, 
contrary to generic situations when the hot parcels as- 
cend and the cold ones descend. Negative Nu may ap- 
pear contradictory, but it could be understood in terms 
of the heat flux through the cross-section at height y, 
Nu(y) = df(y)/dy + J x v(x,y)T(x,y)dx; here J x is the 
integral over the line at height y. We compute this quan- 
tity for various horizontal cross-sections for the flow con- 
figurations corresponding to Fig. [3] (a-f). As illustrated 
in Fig. 0IB), Nu(y) fluctuates significantly in the bulk 
for flow configurations (c,d,e) during the reversals, with 
strong positive values for (d), to strong negative values 
for (e). Note however that Nu(y) is positive near the top 
and bottom plates for all cases. The large fluctuations in 
Nu(j/) in the bulk is a result of the rotation of the central 
role formed after reconnection. We remark that the neg- 
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FIG. 4. (A) Time series of the global Nusselt number. The six 
snapshots of Fig. [3] are marked as (a-f) in this plot. (B) Plot 
of Nu(j/) (the normalized heat transport for the cross-section 
at height y) vs. y during the reversal. The six curves represent 
the six snapshots of Fig. The flow reconfigurations during 
the reversal lead to large fluctuations in Nu(j/) in the bulk, 
ranging from strong positive values for (d) to strong negative 
values for (e). Note however that the Nu(t/) is positive near 
both the plates (y — ¥ 0, 1). 



ative Nu for the two-dimensional convective flow is due 
to strong geometrical constraints faced by the rolls dur- 
ing the reversals, and may not be present in cylindrical 
convection. 

In summary, we show that flow reversals are caused by 
vortex reconnections of two attracting rolls. The flow re- 
configurations during the reversals are due to the nonlin- 
ear interactions among the large-scale modes (1,1), (2,2), 
(1,3), and (3,1). We find large fluctuations in heat trans- 
port during the reversals, which are due to the above 
restructuring of the flow. The (2,2) mode, critical for the 
dynamics of flow reversals, is born after Ra = 2 x 10 7 (for 
Pr = 1), and its strength relative to the (1,1) mode de- 
creases monotonically afterwards. This is the reason why 
flow reversals are observed only for a range of parameter 
values. 

The role of the large-scale modes described in the 
present letter is analogous to the "cessation-led reversal" 
observed in turbulent convection in cylinder (llj as well 
as in dynamo reversals 0], where the quadrupolar mode 
(equivalent to (2,2) mode) dominates the dipolar mode 
(equivalent to (1,1) mode) during the reversal. Future 



work for different geometries, and Prandtl and Rayleigh 
numbers would provide valuable insights that will help us 
build a comprehensive theory of reversals in convection 
and dynamo. 
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